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Abstract 

In this talk, we present our recent next-to-leading order (NLO) nuclear parton distribution 
functions (nPDFs), which we call EPS09. As an extension to earlier NLO analyses, we sup- 
plement the deep inelastic scattering and Drell-Yan dilepton data by inclusive midrapidity pion 
measurements from RHIC in order to reduce the otherwize large freedom in the nuclear gluon 
densities. Our Hessian-type error analysis leading to a collection of nPDF error sets, is the first 
of its kind among the nPDF analyses. 



1. Introduction 

The global analyses of the free nucleon parton distribution functions (PDFs) are grounded 
on the asymptotic freedom of QCD, factorization and parton evolution. These features allow to 
express the hard-process cross-sections formally as 



where /)s are the scale-dependent PDFs, and frij^h+x denote the perturbatively computable par- 
tonic pieces. The factorization theorem has turned out to work extremely well with increasingly 
many different types of data included in the latest free proton PDF analyses. For bound nu- 
cleons factorization is not as well-established, but it has anyway proven to do a very good job 
id Q S H El 01 m describing the measured nuclear modifications o- bound /o- free in deep inelastic 
scattering (DIS) and Drell-Yan (DY) dilepton production involving nuclear targets. This talk 
summarizes our new NLO analysis nfl. 

2. Analysis Method and Framework 

We define the nuclear PDFs through nuclear modification factors Rf(x, Q£) as follows: 

^(,,e 2 )^f ( x,eV TEQ61M fee 2 ), 

where f CTE( ^ 6AM (x, Q 2 ) refers to a CTEQ set of the free proton PDFs [ 8] in the zero-mass variable 
flavour number scheme. The x and A dependences of Rf{x, Qfys satisfying the sum rules, are 
parametrized considering three modifications: Ry{x, Qfy for both valence quarks, R^(x, Qfy for 
all sea quarks, and R^(x, Qfy for gluons. The nuclear PDFs at Q 2 > are obtained as solutions 
to the DGLAP equations using our own NLO DGLAP solver based on a semi-analytical method 
described e.g. in lliil fioll . 
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All cross sections are computed using the factorization theorem and the initial parametriza- 
tion is adjusted to find the minimum of 



x \{a}) = WifXsU*)). where X 2 N (\a)) = {KM) + Y 

For each data set N, the D, denotes the experimental data value with point-to-point uncertainty 
o-j, and Ti is the theory prediction computed using parameters {a}. The pion data suffers from an 
overall ~ 10% normalization uncertainty <x™ rm , and the normalization factor fy e [1 - (f^™, 1 + 
c™™] i s determined by minimizing \ 2 N - The weight factors amplify the importance of those 
data sets whose content is physically relevant, but whose contribution to x 2 would otherwize be 
too small to be noticed by the automated minimization routine we use. 

Besides finding the central set of parameters S° = {a } that optimally fits the data, propa- 
gating the experimental uncertainties to the PDFs has become an inseparable part of the modern 
PDF fits. The Hessian method ll ill , which we use, trusts on a quadratic approximation 

u ' j u 

The eigenvectors Zk of the Hessian matrix H serve as an uncorrected basis for building the 
PDF error sets Sf. These are obtained by displacing the fit parameters to the positive/negative 
direction along Zk such that^ 2 grows by a certain amount Ax 2 from the minimum AkL Using 
these sets, the upper and lower uncertainty of a quantity X can be written e.g. as 

(AX*) 2 * 2[max/min{x(5 i + )-X(5°),X(5p-X(5 () ),0}] 2 (1) 

k 

where X{Sf) denotes the value of the quantity X computed using the set Sf. Requiring each data 
set to remain close to its 90%-confidence range, we obtain Ax 2 = 50. 
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3. Results and Conclusions 

Now, we briefly summarize the main results from the present analysis, starting with Fig. Q] 
where we show the obtained modifications for Lead — the nucleus relevant for the LHC — at 
two scales. Interestingly, even the rather large uncertainty at small-x gluons becomes notably 
smaller in the scale evolution. This is a prediction that might be testable in the future colliders. 

As the DIS and DY data constitute our main data constraints, we display in Fig. |2] some 
examples of the measured nuclear modifications with respect to Deuterium, 

for different nuclei and compare with the EPS09. The shaded blue bands always denote the 
uncertainty derived from the EPS09 error sets. Their size is comparable to the experimental 
errors, supporting our choice for Ay 2 . The nuclear modification for pion production is defined as 

= 1 d 2 N« A "/dp T dy mm=bias ±d 2 a-i An /dp T dy 

dAu ~ <AUi> d 2 N p //d PT dy ~ cPoffdprdy ' 
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Figure 1: The obtained modifications at Q 2 = 1.69 GeV 2 and at Q 2 = 100 GeV 2 . The black lines indicate the best-fit, 
whereas the dotted green curves denote the individual error sets which combine to the shaded bands like in Eq. (TJ. 




where (N co n) denotes the number of binary nucleon-nucleon collisions and pr,ythe pion's trans- 
verse momentum and rapidity. The comparison with the PHENIX and STAR data plotted in 
Fig. [3] shows that the shape of the spectrum — which in our calculation is a reflection of the 
similar shape in Rq — gets well reproduced by EPSOS0 Figure [3] also presents a comparison 
of the EPS09 gluon modifications R^? with the earlier NLO analyses. The significant scatter of 
the curves highlights the difficulty of extracting the nuclear modifications from the DIS and DY 
data alone. Consequently, also the predictions for pion R&\ n differ significantly as is easily seen 
in Fig. [3] This is actually good news as this type of data, especially with better statistics, may 
eventually discriminate between different proposed gluon modifications. 



The shape is practically independent of the set of contemporary fragmentation functions used. 
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Figure 4: The calculated scale evolution of the ratio F^ n /F9 compared with the NMC data 1 17]. 



Attention should be also paid to the experimentally observed scaling-violations and that the 
DGLAP dynamics is able to reproduce them well. Such effects are most transparent e.g. in 
the small-* structure function ratios versus Q 2 , of which Fig. [4] shows an example. To sum- 
marize, the excellent agreement with the experimental data, x 2 IN « 0.79, and especially the 
correct description of the scale-breaking effects — we argue — is evidence for the applicability 
of collinear factorization in nuclear environment. The best fit and all the 30 NLO nPDF error-sets 
are available for practical use from [18]. Also the leading-order EPS09 sets are provided. 
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